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Abstract. Mark Kac gave one of the first results analyzing random polyno- 
mial zeros. He considered the case of independent standard normal coefficients 
and was able to show that the expected number of real zeros for a degree n 
polynomial is on the order of — logn, as n — > oo. Several years later, Sam- 
bandham considered two cases with some dependence assumed among the 
coefficients. The first case looked at coefficients with an exponentially decay- 
ing covariance function, while the second assumed a constant covariance. He 
showed that the expectation of the number of real zeros for an exponentially 
decaying covariance matches the independent case, while having a constant 
covariance reduces the expected number of zeros in half. In this paper we will 
apply techniques similar to Sambandham's and extend his results to a wider 
class of covariance functions. Under certain restrictions on the spectral den- 
sity, we will show that the order of the expected number of real zeros remains 
the same as in the independent case. 



1. Introduction 

One of the earliest results on the expected number of real zeros of the random 
polynomial given by 

n 

(1.1) P„(x)=^Xfex'= 

k=0 

came from Mark Kac [6 . Kac considered the case when the coefficients are assumed 
to be independent standard normal random variables and was able to show that 
the value of the expected number of zeros is on the order of ^ logn, as n — )■ oo. 
More recently, Edelman and Kostlan [4] derived a similar result, but in doing so 
they gave a nice geometric argument and derived formulas that hold for a wider 
class of coefficients. 

A natural generalization of this problem is to assume some dependence among the 
coefficients. Let X^^Xi, ... be a stationary sequence of normal random variables, 
where the covariance function is given by 

r(A:) =E[XoXfc], r(o) = i. 

Under these assumptions, two important results came from Sambandham. The 
first assumes that T{k) — , where p S (0, i) !8^. In this case, it was shown 
that the expected number of zeros is on the same order as when the coefficients are 
independent. The second result assumes the covariance function is constant; that is, 
V{k) = p for any k, where p G (0, 1) [3 [9]. Here, it was shown that the order of the 



2010 Mathematics Subject Classification. Primary 60H99; Secondary 26C10. 
Key words and phrases. Random polynomials, zeros, dependent coefficients. 
This research was partially supported by NSF grant DMS-0706198. 

1 



2 



J. MATAYOSHI 



expected number of zeros is only ^ logn as n — !• oo, or half as many as before. The 
intuitive reasoning for this is as follows. For the first result, the exponential decay of 
the covariance function is fast enough to "simulate" independent behavior. In other 
words, the dependence between the coefficients disappears at a rate fast enough to 
have no effect on the behavior of the zeros. For the second result, however, the 
constant covariance results in the coefficients tending to have the same sign. This 
causes most of the positive zeros to disappear, leaving only half as many zeros as 
before. 

A further question that could be asked is whether or not this result extends to 
a wider class of covariance functions, where the decay rates are between those con- 
sidered by Sambandham. In what follows, we will show that the same asymptotic 
value holds for the expected number of real zeros when the spectral density (which 
will be introduced below) is positive and continuous. If we assume further that the 
spectral density also has one continuous derivative, we will be able to derive an ex- 
plicit value for the order of the error term. Noting that absolute summability of the 
covariance function guarantees the continuity of the spectral density, we will then 
see that behavior similar to the independent case can be expected for covariance 
functions with a wide range of decay rates. The main result is stated as follows. 

Theorem 1.1. Let Pni'x) be the polynomial given in (jl.ip . where the coefficients 
Xo,Xi, . . . form a stationary sequence of standard normals, with covariance func- 
tion r(fc) and spectral density f {(})). Assume that the spectral density does not 
vanish. Letting N{a,P) he the number of zeros of Pn{x) in the interval (a, it 
follows that 

E[iV(-oo,oo)] ~ -logn, for f eC{[-7T,7r]), 

TT 

E[A^(-cx),cx))] = -logn + 0(loglogn), for f e ([-tt, tt]), 

TT 



as n ^ OO. 

Before we proceed any further, a comment must be made. As mentioned on page 
83 of [5], if we consider the function 

x"Pn (^i) = a;" (^Xo + + • • • + X„^^ 

it can be seen that whenever there is a zero of a;"P„ (i) in (l,oo), there is also a 
zero of Pn{x) in (0, 1). Thus, since the distribution of the zeros of the two functions 
are the same, it is sufficient to only look at the interval from (0,1). A similar 
argument works for the negative real line and allows us to restrict our analysis to 
the interval (—1, 1). By then taking twice the result, we will have a value for the 
total expected number of real zeros. 

Our proof will follow that of Sambandham, with some necessary modifications to 
account for the more general assumptions made on the coefficients. The first step 
is to show that there is a negligible amount of zeros on the intervals (0, 1 — -y^^), 
(1 - i^^,l), (-1 + j^,0) and (-1,-1 + l2El2sii), Following this, we will 
then show that the number of real zeros in the intervals (1 — j^p^, 1 — iHiiSIii) and 



REAL ZEROS OF A RANDOM ALGEBRAIC POLYNOMIAL 



3 



(— 1 + '°^ logK _!__!_•) g^j.g gag}j on the order of 77- logn. One significant difference 
from Sambandham's work is that rather than approximating the number of zeros 
with a speciahy chosen function, we will instead derive the asymptotic behavior 
directly from the Kac-Rice formula using the spectral density formulation of the 
covariance function (see below). The result is that the number of zeros in (—1, 1) 
is on the order of log n and, from the comments above, it follows that the total 
expected number of real zeros is on the order of logn. 



2. Preliminary Results 

We will now discuss some results that are variously attributed to Bochner, Her- 
glotz, and Khinchine (see [H|3]). Any covariance function, r(fc), can be expressed 
as 



J —TT 



where F((p) is real, never-decreasing, and bounded. Furthermore, if F{(l)) is also 
absolutely continuous, we have the formula 

(2.1) r(fc) = r e-^^*/((/.)d(/., 

J — TT 

where /(</>) is called the spectral density of the covariance function. A sufficient con- 
dition for the existence of f{(p) is that r(fc) is absolutely summable. Additionally, 
in this case it will be non-negative, continuous, and of the form 

00 

k— — oo 

Now, recalling that N{a, (3) is the number of zeros of Pn(x) in the interval (a, /3), 
the Kac-Rice formula [5] gives the expected number of zeros on this interval as 



(2.2) 
where 



E[7V(a,/3)] = i / ^i^^dx, 



a(.t) = e[p,?(x)] =^^r(fc-j> 



(2.3) i?(a:) =E[F„(x)P;(x)] =^^r(fc-j>a: 



fc+j 

n 

k=0 j=0 

k+j-2 

k=0 j=0 



c{x) ^ E [{p:Ax)f] = E E - 
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Applying we can rewrite these as 

n n 
k=Oj=Q 

/TT n n 

k=0 3=0 

/TT n n 
k=0 ]=Q 

We are now ready to prove our first lemma. 

Lemma 2.1. For the intervals (-1,-1 + '"^^"g" ), (-1 + -[^,0), (0,1 - j^) 
and (1 — iHliHSU^ ^/jg expected number of zeros is O(loglogrt). 

Proof. Following the method of Sambandham, we will define the function 

n n 

n 71 

(2.5) = e-^''^^'' J2 ^''^y' 

k=0 j=o 

1 — j-n+l g—i{n+l)4i ^ _ ri+lgi(n+l)0 



1 - xe-^'l' 1 - ye''"^ ' 

to assist in our computations. As done on page 85 of [8], plugging into our formula 
for A leads to the expression 



A = J H{x,x)f{(j))d(j} 

1 - xe-^'t' 1 - xe^'t' ■' ^ 

Again following [SI, we also get 
(2.7) 



9iJ(a;,y) 



dy 

^ / 1 — x"'^^e~'^^"'^^ 



fi<t>)d<t> 



1 - xe-^'f' 

{n + l)x"e*("+i)"^ (l - xe"l') - (l - a;"+ie'("+i)'^) (-e*"*) ' 
(1 - xe-"f'f 



f{^)d^, 
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and 
(2.8) 

C- 



dxdy 



f{^)d<t> 



y=x 



/-(n + l)a;"e-*("+i)"^ (l - xe"*"^) - (l - a;"+ie"*("+i)'^) (-e-*"^) ^ 
\^ [l-xe-i't'f J 

f-{n + i)a;»e*("+i)'^ (l - xe"'') - (l - x"+ie*("+i)'^) (-e*"*) \ 
' V (1 - xe^^f J 

For X e (0, 1 — i^^TTj) we have 



A 



^ (1 - xe-"^)(l 



and 



1 



< 



^ (1 - xe^^'t'Yil - xe'-'t'Y 

1 r 1 



(1 - a:)2 (1 - a;e-*'^)(l - aie^"^) 



If we consider the quotient '^"^'^ , then 

A 

Plugging into p. 21) . it follows that 



< ^ < 



A) - 1-x 



(2.9) 



E 



0, 1 



loe 



O (log log n) . 



To handle the interval (—1 + j^It^jO) we will substitute in —x, where 



/(0)d0 



A - / H{-x,~x)f{c^)dc^ 

— TT 



/' 



1 + xe-*"^ 
1 



1 + xe"^ 



_^ (1 + xe-*'^)(l + a:e'"^) 



/(0)d0. 



Likewise, 



TT (1 + xe-"^)^ (1 + xe"^)^ 



/(0)d0 



< 



(1 - x)2 (1 + xe-*-^) (1 + xe}'^) 



fi4>)dcf>. 
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Once more we have the inequahty 



VAC - 52 (C 

— ^<U 



1/2 



< 



Applying (|2.2p . we then have 



(2.10) 



E 



-1 + 



loe 71 



O (log log n) . 



Finally, we will consider the intervals (-1,-1 + i^sMli) and (1 - iHsipii^i). 
Similar to inequality derived on page 87 in 8 , for a; G (1 — '°s " ^ i) we have 



VAC - 52 /C 

<'a 



1/2 



< 



Tfc^o Ej^or(fc- j>' 
TLoE-=or(fc-jV^ 



fe+i-2 ~ 



1/2 



< cn. 



Thus, 
(2.11) 



E 



0(log logn). 



For (—1, —1 + iHSipii) -^e will substitute in —x. Notice that since /(0) is con- 
tinuous and non-zero on [— tt, tt], we can bound it from below by a constant ^ > 0. 
We then have 



A > 



2tt 



^ 1 - (_a;)"+ie-*("+i)'^ 1 _ (-x)"+ie*("+i 



1 + xe-^'l' 



1 -I- a^e*"^ 



m ^^(— X 
fe=o 



,2fc 



since /(</>) = ^ in the independent case. Similarly, we can bound f{4>) from above 
by 1^ < 00. This yields the inequality 

M r -{n + l)(-a;)"e-'("+i)"* (l + xe-"l') + (l - (-a;)"+ie-'("+i)'^) e"*"* 



C < 



27r 



(1 + a;e-*'>) 

[n + l)(-a;)"e*("+i)'^ (l -f xe"*') + (l - (-a;)"+ie^("+i)'^) e*"^ 



21 2,^2fc-2 < 



(1 -I- a:e*'^)2 

„2 



fc=0 

It follows that 



Thus, 
(2.12) 



{~xy- 



\2fe 



VAC - 52 /c 



E 



A 



TV -1,-1 



< 



A 



1/2 



< cn. 



log log n 



O (log logn). 



Combining ((2J| - (|2.12p . the result then follows. 



□ 
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3. Proof of Main Theorem 



Now that we have found an upper bound for the expected number of zeros in 
our initial four intervals, we will spend the rest of our time deriving explicit values 
for the intervals (-1 + l2sl£sjl _i + i ) and (1 - 1 - '"«" ). As stated 

^ n ' log n ' V log n ' n ■' 

before, this will be done by deriving asymptotic values for the expressions given in 
(12.21). We will formulate this result as an additional lemma. 



Lemma 3.1. The expected number of zeros for the polynomial Pn{x) in each of 
the intervals (— 1 + ]2E^ML^ — 1 + j^^-^j) and (1 — ^^^^ , 1 — '°s " ^ jg given by the 
following: 

(i) For/eC([-^,7r]), 

log log 11 1 



E 



N 



-1 



logn 



= E 



N 1 



1 _ log log n 



log n 



2ti 



■ los 



(ii) For f ^C^{[-n,iT]), 
log log n 



E 



-1 



logn 



= E 



2ti 



N 1 



1 _ log log n 



\ogn 



log n + O (log log n) 



Proof. Consider the interval (1 - 1 - For a; = 1 - y, define g{y) = 

^ io'°iogn - ^1^0' Af > be chosen such that /(0) < M, for any g [— 7r,7r]. 
Recalling (|2.2p . for A we have 



A: 



1 



(3.1) 



^ (1 - .xe-*"^)(l - xe'"^)' 

„-i(n+l)0 1 i(n+l) 



n+1 



(1 - xe-'-'t'){\ - xe^'f') 



Ai + O 



Our next step is to determine the asymptotics for Ai. We have. 



1 — 2x cos (p + x' 



:fi<P)d<P + 2 



aiv) 



1 — 2x cos (p + X- 



Looking at the first integral yields 

r9{y) 



l-2(l-y) 1-^ + 0(04) 



= 2 



aiv) 



/(</>) 



10 4> + y- 

For / G C([— 7r,7r]) this becomes 

a{y) 



d0 + O(l). 



2/ #1^.0=Marctan^^(^) 
r + 2/^ y 
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while for f ^ ([— 7r,7r]) we have the more detailed representation 

Jo (t> + y 

-'^('^)arctanf^Uo(loglogn). 



y 

For the second integral in Ai, 



1 — 2a; COS + Jj-gj-j^j^i/s 1 — 2a; cos + x 



\giy) 

It follows that 



g{y) y^ + c/)'^ J(^g(^yjy/3 l-2x COS (j) + X^ 

1 



for/eC([-7r,^]), 





^ ^^(0) arctan I 








V y 






^9{y) 






\ y 



o 



Also, notice that 



n J y \giy) 

Combining these results and plugging into p.ip yields 

A ~ MM arctan f for / e C ([-^, tt]), 



(3 2) ^ ^ ^ ^ 

A = MM arctan O f ^ ) , for / e ([-vr, n]). 

y \ y J \9[y). 

Considering B next, 
(3.3) 

gi0 



(l-a;e-^)(l-.e^^)2^^^)'^^ 

+ l)x"e*("+i''^ (1 - xe'"^) 
(l-xe-*<^)(l-a;e^'^)2 '^^'^^'^'^ 

= Bi+B2 + 0{x"+^{\Bi\ + \B2\)). 
Since we will end up showing that the Bi term dominates, we can rewrite this as 

B = Bi+B2 + . 
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To analyze Bi we will split it into two integrals, 

ra<-y^ cos(f)-l + y 



Bi=2 



(l-2(l-2/)cos<^+(l-2/)2) 
cos 4>—l + y 



— of{(t>)d<i>- 



For the first integral we have 

r9(y) 



■■I 



(l-2(l-y) (l-^+om)+{l-yy 



■/(</>)# 



y 



;J{ct>)dct> + 



When / € 7r, tt]) this becomes 
while for f e ([-7r,7r]) 



21'''' ^^(»^#^M,,etanrgMY 

{<P^ + y^f 



r 

= 2 

Jo 

\ y 

For the second integral in Bi, 



„ # + 



(^) 



(where 0o € (0, ^)) 



arctanl^Uo^ ' 



y9{y)J ' 



COS — X 



mdcp 



< 2 



g(^y) (1 - 2xcos<j) + a;2)2 

|cos<^-x|/(0) 



r(g(v 
Jg(v) 



(1 — 2a;cos0 + a;^)^ 



p7T 



~ 2/(0) 

= 
Thus, 



(</,2+y2)2 



# + 2 



[cos /((/)) 

(g(j;))i/3 (1 - 2a;cos(?i + a;2)2 
|cos(/)-x|/(?i) 



del) 



J(a( 



(g(y)y/3 (1 - 2a; cos + a;2)2 



i9{y)r 



Bir^ 


/(O) / 
— arctan 


^9{y) 




y 


K y 


Bl = 


/(o) 

— t;- arctan 


^9{y) 




y2 


V y 



for/eC([-7r,^]), 



arctan + O , for / e {[-.,.])■ 
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Looking at B2 next, 



IB. 



< 



/' 

J — 7r 



(1 - xe-«<*)(l - a;e*'^)2 
(/, + l).r" 



md<t> 



c{n + 1) 



(1 - xe-^'t'){l - xe^'f 

(1 - yr 



y 



where the last hne foUows from our work on A. Now, notice that for y = log " 
we have 

y log n log log n ygiy)' 
Since (1 — y)'^y is a decreasing function as y increases, it follows that for y e 

/ log log n 1 \ 
^ n ' log n ^ 

(3.4) (n + l)(l-y)"2/< ^^(n + l)^i^^ = 0^ ^ 



s'(y) y Kvaiy) 

This also implies that 

a;"+i(Si + B2) < c{n + ly ^ . 

Thus, 

S ~ ^ arctan ( ^) , for / e C ([-tt, ttI), 
(3 5) ^ ^ ^ ^ 

= arctan (^).O(^), for / . C^h. .])• 

Turning now to C, 

(3.6) 

1 r (n + l)^x^" 

r -(n+ l).T"e-*"'^ (1 - a;»+ie*("+i)<^) 
+ ' L (l-.e-4l-^e-^)- ^^^^^'^ 

/TT _ / -i(Ti+l)<;i I i(n+l)(;i\ i „n+l 
(1 - a;e-'"^)2(l - a;e»<^)2 ■' ^ 
= (7i+C2 + C3 + 0(a;"+i(7i). 

Splitting up Ci gives us 

/■s(y) 1 1 
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For the first term we have 
rgiv) 



t-gvu) 

2 / 7 : ^ tt/W^ 



(l-2{l-y)(l-^+0{ct>^))+{l-yf) 
f{4>) 



When / e C([-7r,7r]) it follows that 



/^^^-a^d^^^arctanf^ 



while for /gC1( [-7r,7r]) 

Jo (0"^ + y"^) W / 

r V y / \y^9{y), 



For the second integral in Ci we have 



ngyv))-" f(s) r fU) 

= 2 / 2^<^ + 2 / — 2' 

Jgiv) (1 — 2xcos0 + a;2) J(g(j/))i/3 (1 — 2x cos (/> + a;^) 

/■(sfe))''' 1 r 



= o 

Thus, 
(3.7) 
For Ca, 



'g(y) {4>^+y'^) J(g(y)Y/^ {I- 2x cos (j) + x^) 

1 



{9i.y)) 



Ci ~ ^ arctan for / e C([-7r,7r]), 

2/^ V 2/ / 

= M arctan + O , for / e ([-tt,.]). 



{n + 1) V" , ,,2(1- y?" 

1 - 

while for C3 we have 



^-=£ (i-x:-V)'(:-.e..) j(^w~°("+')" 



iai< 



(1 - xe-*'^)(l - a;e«<^')2 
< c(CiC2)^/2 (by Cauchy-Schwarz) 
{n+l){l-yY 

Also, note that 



d(t) 



X 



y3 y2 
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From p.4p we know that {n + 1)(1 — y)" = O {l/g{y)), which imphes 
(3.8) 



(n + 1)2(1 -jy) 



2n 



1 



2/5^(2/) 



Thus, we can conclude that 



(3.9) 



C^^arctan 



gjy) 
y 

ajy) 
y 



for/eC(h^,7r]), 



(310) 



Combinmg p.2p . p.Sp . and p.9p . for f £ C ([— tt, tt]) we then have the expression 

VAC - 1 
A 2^' 



while for / G ([— 7r,7r]) this becomes 



(3.11) 



VAC - 52 
A 



Pio) 



arctan 



O 



2/(0) 



arctan 



o 



2?/ \g{y) 



Plugging these into (|2.2p gives us 

1 



(3.12) 



E 



N 1 



logn 



,1 



log log n 



27r 



logn, 



for / e C([-7r,7r]), and 



(3.13) 



E 



N 1 



1 _^ log log n 



logn 



27r 



log n + O (log log n) 



for/eCM[-^,7r]). 

To handle the interval from (—1 + " , — 1 + ) we will substitntc in —x = 

\ ' n ' log n / 

-1 + y, where x e [1 ~ j^, 1 - i^sl^sii). We then have 



(3.14) 

A = 

B = 



1 + xe-*"^ 
1 + a;e-*<^ 



1 + xe*"^ 



/('/')d'/', 



-(n + l)(-x)"e*("+i)'^ (1 + a;e*'^) + (l - (-x) 

(1 + le'"''' 



— ( pi{n+l 



)0) e^0\ 



/('^)d<^, 



C = 



-{n + l)(-a;)"e~*("+i)"^ (l + xe"*"^) + (l - (-a;)"+ie-'("+i)'^) e"*"^" 



\ {1 + xe-^'f') 

-{n + l)(-a;)"e*("+i)'^ (l + xe''>') + (l - (-x)"+ie''("+i)'^) e''''\ 
(l + xe^'t'f J 



f{<f>)d4>. 
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Using the fact that /(vr) = /(— vr), we can apply methods shuilar to those from the 
positive case to derive the formulas 



(3.15) 



for / e C([— tTjTt]), and 



A. 
B . 
C ' 



2/(vr) 
V 

y2 



arctan 



arctan 



arctan 



(3.16) 



. 2/(7r) 
A = arctan 

C = ^^^^ arctan 



9{v) 

y 

ajy) 
y 

9{y) 



gjy) 
y 

y 

gjy) 
y 

+ o 

+ 

o 



1 



9{y) 
1 



ygiy)) ' 
1 



y'^g{y)J ' 



for f E ([— 7r,7r]). Combining p.isp and p.l6|) . we then have the same expres- 
sions as before. That is, p.lSp gives 



(3.17) 

while from p. 161) we again have 



VAC - _ 1 



(3.18) 



^/AC - B^ 
A 



arctan 



2/(^) 



arctan 



2M +0(4t 



g{y) 



O 



g{y) 



Plugging these into (|2.2p gives us 

(3.19) E 

for / e C([— 7r,7r]), and 

(3.20) E 
for/eCi(h^,7r]) 



logn 



log log n 1 

N -1+ ^ ^ ,-1 



log n 



— log n + O (log log n) 
2tt 



□ 



We are now able to prove Theorem ll.il 

Proof of Theorem \1.1[ Combining the results of Lemmas 12.11 and 13.11 we have 

1 



(3.21) 



E[iV(-l,l)]^ -logn, for/eC([-7r,^]), 
E[iV(-l,l)] = -Iogn + O(loglogn), for / e C ([-tt, tt]). 

TT 



RecaUing from Section 1 that E[7V(— cx), oo)] = 2E[iV(— 1, 1)], the result then follows. 

□ 
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